EPJ manuscript No. 

(will be inserted by the editor) 



Numerical solutions of thin film equations for polymer flows 

Thomas Salez^'^, Joshua D. McGraw^, Sara L. Cormier^, Oliver Baumchcn^, Kari Dalnoki-Veress^ , and Elic Raphael"'^ 
^ Laboratoire de Physico-Chimie Theorique, UMR CNRS Gulliver 7083, ESPCI, Paris, France 

^ Department of Physics & Astronomy and the Brockhouse Institute for Materials Research, McMaster University, Hamilton, 
Canada 

November 1, 2012 

Abstract We report on the numerical implementation of thin film equations that describe the capillary- 
driven evolution of viscous films, in two-dimensional configurations. After recalling the general forms and 
features of these equations, we focus on two particular cases inspired by experiments: the leveling of a step 
at the free surface of a polymer film, and the leveling of a polymer droplet over an identical film. In each 
case, we first discuss the long-term self-similar regime reached by the numerical solution before comparing 
it to the experimental profile. The agreement between theory and experiment is excellent, thus providing 
a versatile probe for nanorheology of viscous liquids in thin film geometries. 
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Even though hydrodynamics is a well established field, 
several limits and assumptions remain strongly debated. 
For instance, at which length scale and why does the the- 
ory break down? Can one generalize the concept of vis- 
cosity at small scales? What are the details of flow at the 
liquid/solid boundary? How does a droplet spread on a 
given wet or liquid substrate, and what is the role of the 
so-called precursor film? 

Polymeric liquids provide ideal systems with which 
to probe these fundamental questions. Molecular confine- 
ment effects can be investigated because of the large size 
of polymers and the ease with which one can reach system 
sizes that are comparable to that of the molecule. For in- 
stance, the enhancement of effective mobility has been pre- 
dicted [1] , observed [2l|3] , and was related to a reduction of 
the effective entanglement density near surfaces [HIS] . The 
relaxation to bulk mobility has also been observed for spin- 
coated polymer films [6l[7l|8]. Moreover, larger molecular 
weights induce a larger slip at the substrate [51fTU] , which 
enables the study of slippery systems where the 'no-slip' 
boundary condition must be modified [11]. Another fun- 
damental issue is the leveling of a polymer droplet onto 
an identical film. Whereas Tanner's regime of a droplet 
wetting a solid substrate through a thin precursor film 
is well understood P^[T^[T^[T5] . the connection with the 
opposite regime where the droplet is thinner than the un- 
derlying film p!51IT7] should be explored further [T5] . 

A common feature of these problems is that the liquids 
are viscous so that Stokes hydrodynamics may be applied 
|19] . In particular, thin liquid films are well described by 
the lubrication theory through the thin film equations [20|, 



^ e-mail: thomas.salez@gmail.com 



[2TJ[52- However, apart from linearization 
i23.i24] , these particular equations have not yet been solved 
analytically due to their high orders and non-linearities 
|25) . Mathematical advances can nevertheless be found in 
[261127]. 

In the present communication, we focus on a numerical 
approach. In particular, we implement the two-dimensional 
thin film equations for capillary-driven flows using schemes 
inspired by Bertozzi and Zhornitskaya [281129]. After re- 
calling the general forms and features of those thin film 
equations, we focus on two particular numerical situations 
that are compared to experimental systems: the leveling 
of a thin stepped polymer film [501I5T] and the leveling of 
a polymer droplet on an identical polymer film TF. In 
both cases, we present the algorithm, the solution and its 
intermediate asymptotics that we compare to experimen- 
tal data. Those intermediate regimes are of great interest 
for the non-linear thin film equations since they do not 
depend on the precise initial condition. Thus, they should 
be general intermediate solutions ,32] . 



1 General framework 

In this Section, we present the general thin film equations. 
After listing the common assumptions and boundary con- 
ditions, we address the particular cases where the film 
geometry is invariant along one horizontal or angular di- 
mension. 

As in 18, 30 > 31 , we consider polystyrene films above 
their glass transition temperature Tg ~ 100 °C. Extension 
to any other thin viscous liquid is straightforward using 
the relevant parameters. Here, we estimate the following 
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typical parameters [33 1 134 1 135]: height /iq ^ 100 nm, dy- 
namical viscosity ry ~ 1 MPa.s, molecular weight ~ 
100 kg.mol^^, surface tension 7 ^ 30 mN.m~^, density 
p ~ 1 g.cm"'^ and shear modulus G ^ 1 MPa. Let us 
evaluate the typical capillary velocity Vc, Reynolds num- 
ber Re, capillary length 4 and Maxwell viscoelastic time 



Re = 



— ^ 2 /im.min ^ 

hapvc 




(la) 
(lb) 
(Ic) 
(Id) 



With these parameters, and since we observe slow evolu- 
tion of the liquid surface profile over several tens of min- 
utes [T51I501I5T] . we can make the following assumptions: 
we have an incompressible flow of a viscous Newtonian 
fluid where gravity [36,, disjoining pressure ^37i and iner- 
tia are negligible. This flow is well described by the Stokes 
equation: 

VP r]Av , (2) 



combined with the incompressibility condition: 







(3) 



where P and v are the pressure and velocity fields within 
the liquid. In addition, we assume that the lubrication 
approximation is valid, that is: the profile slopes remain 
small in comparison to 1. Finally, we assume 7 and rj to 
be homogeneous and constant. 

As far as the vertical boundary conditions are con- 
cerned, we consider the case of no shear at the liquid-air 
interface: 

d,vU=,, = , (4) 

where z is the vertical coordinate, and we assume a no-slip 
boundary condition at the substrate: 







(5) 



where v^^ is the projection of the velocity that is parallel to 
the liquid-substrate interface. Note that one could easily 
include slip at the substrate by imposing a nonzero wy \z=o 
(see e. g. pTTjl. 

Using the previous assumptions, we derive the thin film 
equations for two particular invariant geometries: trans- 
verse invariance (Case 1) and axisymmetry (Case 2). 

In Case 1, we assume a spatial invariance of the prob- 
lem in one horizontal direction y, which reduces the prob- 
lem to two dimensions. The height of the profile is given 
by h{x, t), where x is the relevant horizontal direction and 
t the time. The pressure is a priori given by P{x,z,t). 
According to the lubrication approximation, we can ne- 
glect the vertical velocities and write: v = v{x,z,t) e^, 
where is the horizontal basis vector. We then project 



and integrate Eq. ([2|, using Eq. ([3|, Eq. Q and Eq. ([5|, 
and find: 

a.p = , (6) 

the pressure P{x,t) is thus invariant in the vertical direc- 
tion, and: 



v{x,z,t) = ^{z^ ^2hz) d^P , 
2?7 



(7) 



which corresponds to the familiar Poiseuille flow. Volume 
conservation requires that: 



dth + dx dz V = 
Jo 



(8) 



Finally, because the pressure does not depend on z, we 
evaluate P at the free surface through the Young-Laplace 
equation. Since the lubrication approximation implies small 
curvatures, the pressure satisfies: 



(9) 



where Pq is the atmospheric pressure. Thus, combining 
Eq. ([7]), Eq. ([8]) and Eq. Q, we get: 



dth+^dx{h''d^h)=0 
3?7 



(10) 



which is the general capillary-driven thin film equation. 

In Case 2, we assume an invariance of the problem by 
rotation about the vertical axis. Therefore, the height of 
the profile is given by h(r, t), where r is the radial coordi- 
nate. The pressure is a priori given by P{r, z, t). According 
to the lubrication approximation, we can neglect the ver- 
tical velocities and write: v = v{r, z, t) e^, where is the 
radial basis vector in cylindrical coordinates. By proceed- 
ing similarly to Case 1, and including the two principal 
curvatures, we obtain: 



dth^ 



irj r 



dr 



-d^h — ^drh 



= 0. (11) 



It is worth stressing that although this axisymmetric ge- 
ometry has been studied in the past through a similar 
equation '38','39J , Eq. (11 1 is more general since it includes 
the two principal curvatures with no far field approxima- 
tion. 



2 Leveling of a stepped polymer film 

In this Section, we focus on the experimental situation 
described in |30|,l31j : a polystyrene stepped film with ini- 
tial heights hi and /12 (see Fig. [ij , which levels above the 
glass transition temperature due to the capillary-driven 
viscous flow. In the following, we present the dimension- 
less mathematical model before describing the algorithm. 
We then characterize the numerical solution. In particular, 
we study the long-term self-similarity of the evolution and 
the sensitivity of this regime to initial conditions. Finally, 
we compare the results to experimental data. 
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(a) 



(b) 




Figure 1. Schematic of the experiment. The ensemble is placed 
on a solid silicon (Si) substrate, (a) As-prepared stepped film 
in the glassy state iSOWSlf . The heights are of the order of hi ~ 
/i2 ~ 100 nm. (b) Above the glass transition temperature, the 
profile h{x, t) levels with time due to capillary- driven viscous 
flow. 



2.1 Mathematical model 

Since this problem is spatially invariant along one horizon- 
tal direction, we refer to Case 1 introduced in Section [T] 
We are interested in the total height of the free surface 
h{x,t) (see Fig. [T]). Let us introduce the natural dimen- 
sionless variables: 



H = 

X = 
T = 



It 



(12a) 
(12b) 
(12c) 



37? /lo 

and the dimcnsionlcss parameters: 

In 

ho 
ho 

where ho is a chosen length scale of reference. Then, we 
non-dimensionalize Eq. ( [lO| ) to obtain: 

OtH + dx (H^d^H) ^0 . (14) 



Hi 

H2 



(13a) 
(13b) 



In addition, we consider the following dimcnsionlcss initial 
condition: 



HiX,0) = 




X <0 
X >0 

X = 



(15) 



P ctra TTi o^" ay 
L cLi cLilifc; LtJi 




AT 


1.0 X 10"^ 


AX 


1.0 X 10"^ 


N 


10000000 


M 


501 


Hi 


1.0 


H2 


1.0 



Table 1. Optimized parameters for the numerical resolution 
of Eq. (fTsll. 



volume conservation, it has been shown in [2^ that this 
method ensures positivity of the height profile H{X, T). In 
the following, we describe the discretization as well as the 
initial and boundary conditions of our algorithm, before 
presenting the integration routine. 

We discretize space through the definition: 



M + 1 



AX,T 



(16) 



for i G [l,Af], where AX is the dimensionless spatial in- 
crement. We introduce the growth rates: 



'H-x,i — 



AX 



T^xxxA — 



AX 

T^Xx.i ~ T^Xx,i-i 



AX 



(17a) 
(17b) 

(17c) 
(17d) 



Then, we consider a continuous-time discrete-space ap- 
proximation of Eq. (14): 



^T'H^ 



AM 



i'^XXXA 



— Ai+iH 



XXX,i+l 



AX 



(18) 



with i g [l,Af — 1], and where we define the auxiliary 
function Ai as: 



Ai 



■W 2 112 



(19) 



for i g [1, M]. Note that for small slopes and small spatial 
increments: Ai ~ , as expected from standard finite 
difference method applied to Eq. ( 14 ) . 



according to Fig.[lJ Eq. ^ and Eq. ([13 1. 
2.2 Algorithm 



The numerical procedure we use is a finite difference method 
developed in [2511^ . It ensures capillary energy and en- 
tropy dissipation as required from [251 . In addition to 



For the initial condition {T — Oj/we use the discretized 
version of Eq. ( 15 ): 



Hi, 

Hi + H2, 
^1 + ^: 



1 < i < 
Af -I- 1 



A-f + 1 



<i< M 



(20) 



Af + 1 
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Figure 2. Numerical solution of Eq. 
sionless times. The initial step is chosen such that Hi = H2 ~ 
1. The inset shows a self-similar replottmg of the data accord- 
ing to Eq. 



For integration, since we have a finite numerical win- 
dow, we choose the following horizontal boundary condi- 
tions at T > 0: 



Figure 3. Evolution of the bump and the dip of the numerical 
solution defined in Eq. for an initial profile with Hi — 

H2 = 1. 



Guided by the symmetry of Eq. (14) and Eq^JTsI, 
we look for self-similar solutions of the first kind |32LI41j 
defined by: 



AmHxxxm 



M 



AX 



(21a) 
(21b) 
(21c) 
(21d) 

(21e) 



F{U) = H{X,T) 
U ^ 



ri/4 



(22a) 
(22b) 



Note that the spatial window size MAX must be chosen 
large enough in order to have relevant horizontal boundary 
conditions. 



We thus replot the evolution of Fig. [2] in this new vari- 
able. The result is shown in inset of the same figure. All 
the curves collapse onto a single curve, demonstrating the 
self-similarity of the long-term evolution. However, this 
intermediate asymptotics appears to be reached only at 
large times, as in [42]. During the short-term evolution, 
the oscillations are observed to grow with time before sat- 
urating. In order to quantify this statement, we study the 
hump and the dip of the leveling profile: 



Finally, we solve Eq. ( 18 1 for any T > using a fourth 



order Runge-Kutta routiiie"|3n] . Time is discretized through Y\, (T) 

T = jAT, where j G [1, N], and where AT is the dimen- 
sionless temporal increment. Note that, according to [38 , Yd{T) 
we should have at least AT ^ AX'^ due to the general 



max[i? (X, T)] - ni&x[H{X, 0)] 



min[i7(X, T)] - mm[H{X, 0)] 



(23a) 
(23b) 



orders of Eq. ( 14 ) . The typical numerical parameters after 
optimization are summarized in Table [T] 



2.3 Results 



In this part, we present the numerical solution of Eq. ( 14 1 
using the algorithm presented above. We characterize the 
long-term self-similarity of the evolution as well as the 
robustness of this regime with respect to variations of the 
initial condition. 

As expected, the initial step levels due to the gradients 
in capillary pressure. The numerical solution H{X, T)—Hi 
is plotted at different dimensionless times in Fig. [2] The 
initial step has been chosen such that Hi = H2 — 1- In 
addition, we observe spatially damped oscillations of the 
free surface. 



Defined in this way, the bump and the dip are the extremal 
heights above and below the as-prepared step profile, re- 
spectively. The temporal evolutions of lb and Id for a film 
with initial condition Hi ~ H2 — 1 are shown in Fig. |3] 
They saturate after T ^ 0.0001, indicating that the self- 
similar regime is reached. 

In order to understand the stability of the self-similar 
regime described in the previous Section, we turned the 
sharp initial profile to a smoother Fermi-Dirac function: 



Hx{X,0)^Hi 



Ho 



1 +exp (-f 



(24) 



where A is the width of the transition region. Note that it is 
necessary to adapt the spatial increment so that the condi- 
tion A > AX is fulfilled. The results for A = 0.1 and A = 
are compared in Fig. [4] Both initial conditions converge 
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Figure 4. Numerical solution for a smooth initial condition 
given by Eq. Q with ffi = Ha = 1, A = 0.1 and AX = 0.01. 
The profile is plotted at several dimensionless times, in self- 
similar variables according to Eq. \2t^ . For comparison, the 
solution from a pure Heaviside initial condition (X — 0) has 
been plotted as well. 



to the same self-similar regime of Eq. ( 22 1 . The smaller A 



is, the faster the self-similar regime is reached as expected 
according to the driving Laplace's pressure of Eq. Q. The 
asymptotic regime is thus attractive and robust with re- 
spect to a variation of the initial profile at constant bound- 
ary limits. This result is important for comparison to ex- 
periments since it validates a posteriori the contradictory 
use of a sharp stepped initial condition in Eq. ( 15 1 within 



the lubrication approximation of small slopes underlying 
Eq. (14 1. Moreover, this tells us that experimentally fab- 
ricatea samples with their inherent imperfections should 
still approach the theoretical self-similar regime, as dis- 
cussed below. 



according to Eq. ( 12 ). For the data shown in Fig.[5j we find 
7/77 = 1.5 ± 0.3 jjm.min""'^ at 140°C [3T], which compares 
well with the tabulated values through the WLF 

model [H] . This agreement validates the numerical results 
and the methodology itself. Moreover, the technique ap- 
pears to provide an accurate measurement of the capillary 
velocity of the liquid. 



3 Leveling of a polymer droplet on an 
identical film 

In this Section, we address a second problem inspired by 
experiments; the leveling of a thin liquid polymer droplet 
on an identical film 18J. We assume that the initial shape 
of the droplet is a spherical cap of height dc and radius of 
curvature Tc, and that the underlying film of thickness e is 
infinitely wide (see Fig. |6). Whereas Tanner's regime of a 
droplet wetting a solid substrate is well understood by in- 
voking the presence of an infinitesimal precursor film |121 
[T3 l [T4 j [T5] , the connection towards the opposite regime [TBI 
[T71ITS] . where the droplet is thinner than the film, is of 
interest. In the following, we establish the dimensionless 
mathematical model before presenting our numerical algo- 
rithm. Then, we characterize the resulting numerical solu- 
tion. In particular, we study the long-term self-similarity 
of this evolution. Finally, we compare the results to exper- 
imental data. 



3.1 Mathematical model 

Since this problem is axisymmetric, we refer to Case 2 in- 
troduced in Section[l] We are interested in the total height 
of the free surface hjr, t) (see Fig. [6]) . Let us introduce the 



2.4 Comparison with experiments 

In order to demonstrate the interest in such a numeri- 
cal solution, we now compare it with experiments. Details 
can be found in [501131) . Note that the self-similarity of 
the experimental profiles has been demonstrated in [31 . 
Figure [5] shows a comparison between the numerical solu- 
tion of Eq. (14) and an experimental profile. The sample 
is a polystyrene stepped film {hi — h2 = 89 nm) with 
= 118 kg.mol"^, annealed at 140°C. As we can see, 
the agreement is excellent, with residuals being less than 
one percent of the data range. The vertical scaling pa- 
rameter is determined by according to Eq. ( 12 ) and 



Eq. ( 13 1. Thus, if the capillary velocity 7/77 is unknown at 
the considered temperature, the single fitting parameter 
is the horizontal stretch defined by: 



U = 



X 

ri/4 

3iT_ 



1/4 



^1/4 



(25a) 



(25b) 




Figure 5. Fit of the numerical solution of Eq. ,14-) to experi- 
mental data ( atomic force microscopy) , in self-similar variables 
according to Eq. \2S^ . The system is a polystyrene stepped film, 
with Mm — 118 kg.mol~^ and h\ = h2 — 89 nm, leveling on 
a silicon substrate !S0l\31f . The temperature of the sample is 
maintained at 140° C. 
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Figure 6. Schematic of the experimental configuration in \1S^ - 
The ensemble is placed on a solid silicon (Si) substrate, (a) 
A thin polystyrene droplet, of initial height dc and radius of 
curvature r^, is placed onto an identical film of thickness e, 
in the glassy state, (b) After heating above the glass transition 
temperature, one observes the capillary leveling of the profile 
h{r,t). 



natural dimensionless variables 
R 



H = 
T = 



r 

ho 
h 

ho 

It 
377/10 



and the dimensionless parameters: 

Rc = 
E 



dc_ 
ho 

ho 

e 

ho 



(26a) 
(26b) 
(26c) 

(27a) 
(27b) 
(27c) 



where ho is a chosen length scale of reference. Then, we 
non-dimensionalize Eq. (11) to obtain: 



dr- 



RH 



1 



In addition, we consider the following dimensionless initial 
condition: 



H{R,0) 



E 
E, 



D,-R, + ^R^-R\ 



R<A 
R> A, 

(29) 



according to Fig. [6j Eq. (26) and Eq. (27). We introduced 
the distance A from the center to the initial point where 
the flat film intersects the droplet profile: 



A = ^2R,D, - D'i 



(30) 



3.2 Algorithm 



We use a similar algorithm as the one introduced in Part 2.2 
In the following, we describe the discretization as well as 
the initial and boundary conditions, before presenting the 
integration routine. 

For spatial discretization, we define: 



M + 1 



AR,T 



(31) 



for i € [1,M], where AR is the dimensionless spatial in- 
crement. We introduce the growth rates: 



^ _ ^i+i — 'Hi 



n 



RR.i 



Tin 



RRR,i 



AR 

'^RR,i ~ 'HrRA-1 

AR 



(32a) 
(32b) 

(32c) 
(32d) 



Then, using Eq. ( 19 ), we consider a continuous-time discrete- 



space approximation of Eq. ( 28 ) : 



+ 



■^i^RRR,i ~ ■^i+l'H-RRR.i+l 

AR 

■^i^RR,i " ■^i+iT~l-RR,i+l 
(7^^3)AR2 

■^i+lT~l-R,i+l ^ ■^i^R,i 

{i - iYAR^ 

■^iT~l-RRR,i _ -^i^R,! 

{i^i)AR {i ~ 3)3 AR^ ' 



(33) 



with i€[4,M - 1]. 

At T = 0, we use the discretized version of Eq. ( 29 ) as 
an initial condition: 



(28) n, = 



'e + Dc-Rc+ a/-Rc - - 3)2Z\i?2, 



E, 



i < k 
i > k , 
(34) 



with: 



fc = 3 



AR 



(35) 



where we introduced the floor notation [ J . 

For integration, since we have a finite numerical win- 
dow, we choose the following horizontal boundary condi- 
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Parameter 


Value 


AT 


1.0 X 10"^° 


AR 


1.0 X 10"^ 


N 


100000000 


M 


300 


Dc 


0.1 


Rc 


2.0 


E 


1.0 



Table 2. Optimized parameters for the numerical resolution 
of Eq. pSl). 



tions at T > 0; 

-Hi 

n2 

T~Lrrra 
T~Lrr.m 



dr-Hs = 2 



A5 (36a) 

H5 (36b) 

Hi (36c) 

-^fl,5 (36d) 

-'^rrr,5 (36e) 

(36f) 

(36g) 

AR • 



The last equation comes from the parity of the profile in 
i' = i~ 3 and the identity of the two principal curvatures 
at i' = 0. The value i' = corresponds to i? = 0, since 
we need two additional points to define the central growth 
rates. Again, we note that the spatial window size MAR 
must be chosen large enough in order to have relevant 
horizontal boundary conditions. 

Finally, we solve Eq. ( [33| for any T > using a fourth 
order Runge-Kutta routine 40 . Wc discretize time through 
T = j^T, where j £ [i,N] and where AT is the dimen- 
sionless temporal increment. Note that, according to |38 , 
we should have at least AT ^ AR"^ due to the general 



orders of Eq. ( 28 ) . The typical numerical parameters after 
optimization are summarized in Table [2] 



3.3 Results 



In this Part, we present the numerical solution of Eq. ( 28 1 
using the algorithm introduced above. Then, we charac- 
terize the long-term self-similarity of the evolution and the 
convergence to this regime. 

As expected the initial drop spreads due to the large 
gradients in capillary pressure. The numerical solution 
H{R, T) is plotted in Fig. [t] at different dimensionless 
times. Here, the initial conaition is a spherical cap of 
dimensionless height = 0.1 and radius of curvature 
i?c = 2, on a film with height E = 1. Note the presence of 
oscillations of the free surface, and especially a dip, as in 
Fig. [2} 



7=0.00000001 

T=0.0001 

T=0.001 

T=0.005 

T=0.01 

T=0.05 




Figure 7. Numerical solution of Eq. (28), after several di- 
mensionless times. The initial condition is a spherical cap of 
dimensionless height Dc — 0.1 and radius of curvature Rc = 2, 
on a film with height E = 1. 




w 



'28), after 



Figure 8. Long-term numerical solution of Eq. 
several dimensionless times, in the self-similar variables of 
Eq. (37| j. The initial condition is a spherical cap of dimen- 
sionless height Dc = 0.1 and radius of curvature Rc = 2, on a 
film with height E = 1. 



Guided by the symmetry of Eq. (28) and Eq. ( 29[ ) 
we look for self-similar solutions of the second kind ^32, 
defined by: 



G{W) = 
W = 



H{R,T)-E 
H{0,T) ~ E 
R 

[H{0,T)3 T]i/4 



(37a) 
(37b) 



We thus plot the long-term evolution of the droplet of 
Fig. JT] with these new variables. The result is shown in 
Fig. |8j All the curves collapse onto a single curve, demon- 
strating the self-similarity of the long-term evolution. In 
order to characterize the latter, we introduce the addi- 
tional dimensionless volume V with respect to the volume 
of the underlying film. It is constant by incompressibility 
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1x10 



0.0001 



Figure 9. Crossover between the two regimes of a levehng vis- 
cous droplet onto an identical film. Data was computed accord- 
ing to the numerical solution of Eq. 
Ho introduced in Eq. i3!X. 



'28) and the function of 




• Experiment 
— Theory 



Figure 10. Fit of the numerical solution of Eq. .28) to 
the experimental profile (atomic force microscopy), in self- 
similar variables according to Eq. [g?] ). The system is a thin 
polystyrene droplet with M^, = 118 kg.mol~^ leveling on an 
identical film, at 140 ° C '181. The dimensionless parameters 
are: T = 105, Dc = 0.98, Rc = 67.5 and E = 1. 



of the liquid, and is defined by: 

/>oo 

V = 27r / dRR [H{R, T) - E] , (38) 
Jo 

which can be evaluated at T = since we know the ini- 



tial profile. Using Eq. (371 and Eq. (38), we derive the 



temporal evolution of the droplet heig 

3 / rr \ 2 



Ho 
E 



1 - 



Ho 
E 



It Ho = H{0,T) 
V2 



where: 



a = 27r / dW W G{W) 
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(39) 



(40) 



The last quantity is a geometrical factor depending only 
on the dimensionless height and radius of curvature 
Rc of the initial spherical cap. However, Eq. (39) and thus 



self-similarity of Eq. ( 37 ) appear to be only true at large 
times, as in During the short-term evolution, the drop 
height does not change significantly since the Laplace pres- 
sure gradients are much greater near the initial contact 
line. Therefore, we expect a crossover between a transient 
regime of constant height during which the contact line 
at the periphery of the droplet becomes smooth - and the 
influence of the initial condition is vanishing - and the in- 
termediate asymptotics of Eq. ( 37 ) . The full evolution of 
Hq is shown in Fig. [9] It confirms the existence of two 
regimes. After T ~ O.I, the evolution is self-similar. More- 
over, the long-term slope of the curve in Fig. [9] is self- 
consistent with the a and V values of Eq. ([39|, for each 
configuration that we tested numerically. 



3.4 Comparison with experiments 

In order to validate the previous computational results, 
we compare the numerical solution to experimental pro- 
files. The system is a thin polystyrene droplet with = 



118 kg.mol""'^ leveling on an identical thin film, at 140 °C. 
Details can be found in |I8j , where we use the present nu- 
merical solution to establish the connection between Tan- 
ner's regime of a droplet spreading onto an infinitesimal 
precursor film [T2lI3, 14, 15J and the opposite regime of a 
small droplet leveling into a thicker film. A typical result 
is shown in Fig. 
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Note that the self-similarity of the 
experimental profile has been verified. There is no free pa- 
rameter in this comparison since the capillary velocity has 
been extracted from the leveling profile of a stepped film 
(see Section [2]) with the same liquid and temperature [?T] . 



Conclusion 

We presented the numerical implementation of capillary- 
driven thin film equations that describe the evolution of 
viscous films in two-dimensional configurations. In both 
cases, we presented the algorithm, the solution and its in- 
termediate asymptotics that we compared to experimental 
data. 

The first application was the leveling of a step at the 
free surface of a thin polymer film. In this case, we demon- 
strated the self-similarity of the long-term evolution and 
characterized the robustness of this intermediate asymp- 
totics of the first kind. Then, we compared the numerical 
solution to experimental profiles. The excellent agreement 
validates this technique as a useful viscometer for polymer 
melts. 

The second application that we discussed was the lev- 
eling of a liquid polymer droplet onto an identical film. 
In that case, we found another type of self-similar evolu- 
tion. This intermediate asymptotics of the second kind, 
as well as the convergence to this regime, were character- 
ized as in the case of the leveling of a stepped film. The 
evolution of a given droplet appeared to crossover from 
a transient regime of constant height where the initial 
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contact line becomes smooth to the long-term self-similar 
regime. The agreement with experiments is also excellent. 
Whereas Tanner's regime of a droplet on an infinitesimal 
precursor film is well understood [l2 l [T3 l [T4 l [T5j . the nu- 
merical approach presented here offers a way to explore 
further the regime where the thickness of the underlying 
film is no longer negligible with respect to the height of 
the droplet [T5] . 

In near future, the computations may be extended 
to other geometries, fluids, and length scales. Disjoining 
pressure at the smallest scales [37] or gravity at much 
larger scales [36] may also be included. It is our hope that 
the simple approach presented here may be more gener- 
ally used to probe rheology and boundary conditions on 
the nanoscale, and facilitate an understanding of hydro- 
dynamics on length scales where continuum treatments 
break down. 
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